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Abstract. - Many models for chaotic systems consist of joining two integrable systems 
with incompatible constants of motion. The quantum counterparts of such models have a 
propagator which factorizes into two integrable parts. Each part can be diagonalized. The two 
eigenvector bases are related by an orthogonal (or unitary) transformation. We construct a 
random matrix ensemble that mimics this situation and consists of a product of a diagonal, an 
^—s . orthogonal, another diagonal and the transposed orthogonal matrix. The diagonal phases are 

0\ ' chosen at random and the orthogonal matrix from Haar's measure. We derive asymptotic results 

f^ , (dimension A'' ~+ oo) using Wick contractions. A new approximation for the group integration 

f^ ■ yields the next order in 1/A'^. We obtain a finite correction to the circular orthogonal ensemble, 

f^ ' important in the long-range part of spectral correlations. 



Ft, . 

^1 It is usually assumed that the spectral fluctuation properties of classically chaotic sys- 

• ^ ■ terns coincide with those of the corresponding canonical ensemble of random-matrix theory 
rS I ("quantum chaos conjecture"). There is abundant numerical evidence for the conjecture. In 
i^ ■ addition, several approaches have aimed at an analytical proof for the conjecture, see Section 
5.9 of the review |l|. Naturally, these approaches have addressed generic systems. However, 
many of the commonly studied chaotic systems have a very particular form: They can be 
divided into two parts each of which is integrable, albeit in different coordinates. Such sys- 
tems have been called semi-separable 0]. Typical cases are: The quarter stadium, kicked 
systems such as the kicked rotor, chaotic Jung scattering maps for integrable Hamiltonians 
[^ , or Lombardi's approximation for Rydberg molecules Q] . The spectral fluctuations of such 
semi-separable systems have been found to be essentially consistent with the quantum chaos 
conjecture. In view of the special nature of these systems, that finding is somewhat surprising. 
In the present letter, we construct two random-matrix models (one for the unitary and one for 
the orthogonal case) for semi-separable systems which we can solve analytically. The models 
take account of very specific features of semi-separable systems which are not reproduced by 
the circular ensembles of random-matrix theory. With the help of these models, we show why 
semi-separable systems nearly follow the predictions of standard random-matrix theory, and 
we predict quantitatively the deviations that are typically expected for such systems. In order 
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to perform the ensemble averages, we develop a novel technique of approximate integration 
over the orthogonal or the unitary group in the limit of large matrix dimension. 

Semi-separable systems are computationally relatively simple: One uses a surface of section 
that separates the two integrable parts. For kicked rotors one gets an intertwining of gauge 
transformation and Fourier transformation if (time) sections are taken before and after each 
kick. For scattering maps one obtains transport with the full and the free Hamiltonian if 
the cut is placed in the asymptotic region. In the case of time-independent systems with 
two degrees of freedom, one can often choose a Poincare section such that both parts become 
integrable scattering systems when the other part is replaced by a scattering channel. The 
quantum Poincare map (QPM) is then the product of two Jung scattering maps ||, ||. The 
latter can be unitarised [g|. In particular cases such as the quantum-defect theory of the 
Rydberg molecule, the division may relate to a variable such as angular momentum that 
has a discrete quantum number, leading to a unitary QPM Q. For surfaces of section, the 
relevant ensembles to implement the quantum chaos conjecture are the circular orthogonal 
and the circular unitary ensemble (COE and CUE). 

Typically, the time-evolution operator or the QPM of a semi-separable time-reversal 
invariant system consists of a sequence of four operations: The orthogonal transformation 
O which takes us from the basis where one part of the propagator is diagonal to the one 
where the other part of the propagator is diagonal, followed by a diagonal phase matrix 
Di — diagje^*''"}, this by the inverse transformation matrix O^ and this by another diagonal 
phase matrix D2 = diagje^*^"}. For the propagator T' this yields the N x N matrix T' = 
D20^ DiO. The time-evolution operator can equivalently be written as 

T = D^^^^T'dI^^ = D^'^O^DiODy'^ . (1) 

If time-reversal invariance is broken in either of the two sub-systems, the orthogonal matrix 
O, and its inverse O"^, have to be replaced by a unitary matrices C/, and C/^ — U^^. The 
ensembles are then generated by assuming that the phases in Di and D2 are independent 
random variables, i.e., that they obey Poissonian statistics, and that O and U are members of 
ensembles defined by the Haar measure for the orthogonal and unitary matrices, respectively. 
We denote the resulting two ensembles for T by DODO and DUDU, respectively. Note that 
such and similar composed ensembles of random matrices have been studied numerically in 
Ref.§. 

We determine the spectral properties of DODO and DUDU by first computing the spectral 
form factors K{m) — (jtrT™!^) for both cases. From K{m) we can calculate any 2-point 
statistics. For example, the number variance Yp'[L]N) is given by 

m=l ^ ^ 

Here L is the length of the spectrum over which the variance is calculated. To leading order 
in 1/A^, the asymptotic result for DODO is given by 

S^ODo(i; N) = Y^lo^m N) + \ sin^TTL/N) + O il/N) , (3) 

TT 

and similarly for DUDU. 

In order to compute K{m), we must calculate averages of powers of T. For the quadratic 
term, we have 

(rfc,T;,}DODO = (e*^'"=+'"-''''-'''))^^(e2*«"-«'))e(OfcaOiaOpfc09b)o . (4) 

a.b 
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The averages ()^ and ()j are trivial and yield (e"'*-' "''),, = S^^i and correspondingly for ()j. 
This gives (Tfc,Tp*^)DODO = A(/cZ,pg) J2a{OkaOiaOpaOqa)o- The factor A{kl,pq) equals unity 
if the pairs of indices kl and pq coincide (so that either k — p and I = q ov k = q and I = p) and 
vanishes otherwise. The difficulty lies in calculating the average ()o over the orthogonal group 
0{N). This is done approximately. We begin with the terms of leading order in \/N . We 
assume the matrix elements to be independent variables, and we replace the group integration 
by a Gaussian average ()g over the space of real N x N matrices, 

For the variances, this yields the correct expression {Oi.jOki)G — jfSikSji- In this approxi- 
mation, averages of monomials of matrix elements of O can be calculated by standard Wick 
contraction rules for Gaussian averages. Because of the symmetry of the matrix T, the only 
two non-equivalent non-vanishing quadratic moments of T are 

(Tiir*i)DODO ~ 3/7V, (Ti2r*2)DODO ~ l/N . (6) 

Comparing this to the result for COE we find that only the diagonal term differs to leading 
order: The COE result is 2/N. For K{1), this yields K{1) = 3 which leads directly to 
Eq. (|^). Higher moments of T are not affected by the difference because when compared to 
non-diagonal contractions, the contribution of diagonal contractions is suppressed by a factor 
1/iV. 

In order to obtain the subleading terms in 1/A^, we must improve our approximate group 
integration. We note that the Gaussian measure in Eq. (^) yields orthogonality of the matrix 
0,nn only for quadratic terms, {{0"'"0)ki)G = ^k,i- For higher polynomials in 0„i„, corrections 
in the first subleading order occur. For example, 

{{0'^0)k,i, (O^O)fe.i. • • • {0^0)k„jjG = 6k,i,6k,i, ■ ■ ■ 5k„j„, + O il/N) , m > 2 . (7) 

We observe that as A^ — > oo, our approximate integration in Eq. (H) is increasingly restricted 
to the neighborhood of the group manifold. To obtain the desired improvement, we introduce 
an extra weight factor w{0) which is chosen such that the orthogonality relations (R) become 
exact for larger values of m, or the error terms O (1/-/V) are replaced by terms O (l/A^^j, or 
both. We define the new "average" ()„, 

(/(O)).. = {I{0)w{0))g ■ (8) 

To be useful, the weight function k;(0) has to be of low order in the matrix elements of O 
and should have a limited number of free parameters to adjust. We use the ansatz 

w{0)^A + BiY{0'^0) + Ctr{0'^00'^0)^D[tr{0'^0)f . (9) 

This expression is of fourth order in O and contains all essentially different terms up to this 
order. The use of traces of powers of O^O guarantees the symmetry under interchange of 
indices. Moreover, it causes w{0) to be invariant under left and/or right multiplications with 
an arbitrary orthogonal matrix O', 

w{0) = w{00') = w{0'0)., O' G 0(iV) . (10) 

We note that O is not a member of the orthogonal group. The invariance property (|lO| ) is 
desirable because it guarantees that the measure on the group manifold remains unchanged 
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except for a constant, and that areas close to this manifold are not affected drastically because 
they will be mapped onto other areas near the manifold. We determine the free parameters 
A, B, C, D by imposing the correct orthogonality properties for terms in O^O of order 
(normalization), 1, and 2, 



(1)^ = 1 , (11) 

((O^O),,)^ = ,5,,, , (12) 

((0^0),,,(0^0)fc,z)^ = S.JkA ■ (13) 

Using the approximate integration, we find for the left-hand side of Eq. ( [l3|) the form 
UA, B, C, D, N)5,^,5kj+fi{A. B, C, D, N){6,^k5jj+5,^^k)- This yields h{A, B, C, D, N) = 
1 and f4{A,B,C,D,N) = 0. Together with Eqs. ( [ll|jl^ ) one finds four linear equations for 
the four unknowns A, B, C, D. These were solved using Mathematica 4. We obtain |(^ 

A^l , B = —, C= -^ -, D^^-^ ^. (14) 

4 ' 2 ' A{N^ + N-2) A{N^ + N-2) ^ ' 

The negative signs in A and C imply that w{0) is not positive definite and, thus, is not 
a proper measure. Yet the constant on the group manifold is positive, namely 1 + N^/A. 
Therefore, on the part of the integration area near this manifold where the Gaussian factor is 
not small, the weight w{0) is also positive. This fact reflects the invariance condition (pO|). 

In addition, we get the bonus that with the new integration (^ replacing the Gaussian 
integration (||) the correction terms on the right-hand side of Eq. 1^) for to > 3 are O [^). 
We have checked this statement explicitly for m = 3,4,5, and we conjecture that it holds 
for arbitrary to > 3. We have also checked that averages ()u, over monomials containing 
the components Oia of a single vector (fixed a) up to 10th power are correct to subleading 
order 0{^). We also applied our procedure up to products of three indcntitics {0^0)ij in 
the postulated relations (O) - ([l3h with a unique solution for 6th order polynomial ansatz for 
w{0). For 7TJ > 4, no improvement over the 4th order case (p|,p^ was found in the approximate 
orthogonality relations Eq. (M) (again within 0{^)) for terms of power higher than 6; only 
the latter became exact. 

We calculate the form factor for DODO by averaging over all the Poissonian phases, and 
by replacing the group integral by ()^. For ttt, = 1, 2, this yields 

OAT 

^DODO(I) = E(Ofep)» - Jfp^ , (15) 

i^DODo(2) = 4 Y. {ol,,ol^p,ol^pAp.U-^T.^ol^A.)-' + T.<ol;)^{i&) 

ki,k2,pi,P2 k,pi,p2 k,p 

= 4-l + 0(iV-2) . (17) 

The first moment is exact as the averages involved are of the type to which the weight factor 
was adjusted. For the second moment we stipulate that Qw correctly reproduces to subleading 
order, not only the averaged expressions mentioned above, but also the first term on the r.h.s. 
ofEq. (0). 

(^)The routines for averaging of symbolic sums of monomials of matrix elements over orthogonal and unitary 
groups using the approach described here can be obtained from one of the authors by e-mail request on 
prosenOf iz.uni-lj . si. 
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Fig. 1 - Deviation of the number variance of DODO(A'^ — 50) from COE{N — 50) as determined 
by numerical Monte Carlo simulation (points), compared to the theoretical prediction with 1/N 
correction (nsl) (solid curve) and without 1/A'^ correction (A'^ = oo) (dashed curve). 



We determine the number variance from Eq. (Q), terminating the sum after the second 
moment. We expect the result to give a good account of the long-range correlations. We can 
calculate the COE result for finite N in the same approximation either in the conventional 
way or using similar techniques, see below. For the difference we obtain 
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(18) 



To lowest order in 1/N, this formula agrees with Eq. (0) and is exact, because only the first 
moment of T has a correction of order zero in 1/A^. We expect that the terms of next order 
in 1/A'^ yield a good representation of the long-range behavior in L of the number variance. 
Finite-size effects for the short-range behavior are difficult to predict, except that we expect 
them to be small on general grounds, at least as long as the range L is short compared to the 
matrix dimension N. 

We used numerical calculations to test our results and to look for finite-size corrections 
at short distances in the spectrum. Fig. H shows the difference of the number variances for 
DODO and COE for three cases: In the limit N ^ oo and for A^ = 50 (in both cases, DODO 
results were obtained from Eq. dlq), and the numerical result for A^ = 50. We find excellent 
agreement (even at short distances) confirming our intuitive argument. Since the number 
variance is at short distances not very sensitive to the changes considered in this paper, we 
show in Fig. g the difference between the integrated nearest-neighbour spacing distributions 
I{S) = J^ dsP{s) for DODO and COE for A^ = 25, 50, and 100. For purposes of reference, 
we have also included a comparison with the Wigner surmise. This shows how small the 
correction really is, which for A^ = 100 disappears in the statistical noise of 2 x 10^ spacings. 
We mentioned that the ensemble we use has been studied numerically before Q . The authors 
only looked at fairly short spectral ranges. This is why they failed to observe the zero order 
deviation and the A^ -dependence. 
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Fig. 2 - Deviation of the integrated level spacing distribution of DODO(7V) from COE(7V), 
-fDODO(iv)('S')— /coE(]V)('S'), as determined by numerical simulations for A'^ — 25, 50, 100 (thin, medium, 
thick curve, respectively). The deviation of the Wigner surmise from COE(oo) is shown by the dashed 
curve. 



For the unitary case, the approximate group integration can be developed in fuU analogy 
to the orthogonal case. For the weight function Wc{U) (where [/ is a complex A^-dimensional 
matrix) we use the ansatz Wc{U) ^ A^ + B^MU^U) + CctT{WUUW) + Dc[tT:{WU)]^ and 
find 

(19) 



1 



B,^N 



Cr = 



2{N^~1) 



Dr = 



2(iV2 - 1) 



Again, the integral of any monomial {f{U))u is aproximated up to the first two (leading and 
subleading) orders in 1/N by applying complex pairwise contraction rules to the polynomial 
f{U)wc{U). The results for the first two moments for DUDU are 



K, 



DUDU 



(1) = 



2N 
N+1 



i^DUDu(2) = 2-l + o(^ 



(20) 



We have similarly calculated the leading and subleading terms for the COE and the CUE. 
The relevant matrices are represented as U^U and U, respectively. The results coincide with 
the leading terms of the 1/A^ expansion of the exact expressions. 

Fig. 1 shows that the DODO spectrum is less stiff than that of the COE. A corresponding 
statement applies to DUDU versus CUE. The differences are caused by the Poissonian statis- 
tics of the phases rjk and ^a- To see this, we replace either of the two diagonal matrices Di 
and D2 in Eq. (nf) by another diagonal matrix containing the eigenvalues of the COE (or of 
the CUE). It is straightforward to show that the resulting ensemble is identical to the COE 
(the CUE, respectively). In the opposite limit, we keep D2 (Poissonian) and replace Di by 
diag{exp(ian)} (picket fence model). We denote the associated first moment by -^dodo(-'^)- 
Averaging over ^„, summing over the n's, and calculating the averages of the O's with the 
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help of the technique described above, we find 

N N N 

DODOW A^ A^ f\ \ )) A^\ Ira /n/u ^ j^ ^ (iV + 2) sin^ a/2 ^ ' 

m—1 n—1 1=1 ^ ^ \ / / 

As we deal with a fourth-order monomial, this is an exact result. Eq. ( |2l| ) yields an interpo- 
lation for the long-range level statistics from the Poissonian value K{1) = N (this applies in 
the limit of a degenerate spectrum where a = 0) to the COE value Kcoe{^) — 2N/{N + 2) 
(with error 0{1/N)) when a becomes bigger than 2-k/N. This result suggests that DODO 
or DUDU provide a worst-case scenario for possible deviations from the canonical random- 
matrix predictions of COE or CUE. 

In summary, we have introduced random-matrix ensembles for semi-separable systems. 
These allow us to understand why semi-separable systems follow the quantum chaos conjecture 
very closely, and to predict quantitatively the deviations. On scales L which are of the order 
of 10% of the spectral range N or less, the spectral fluctuations for each ensemble are very 
similar to those of the corresponding circular ensemble. Beyond this range, the deviations 
reach values of up to 0.2 or so. On a relative scale, these deviations are small (we recall the 
logarithmic increase of S^). On an absolute scale, they are not. Furthermore, they show a 
well-understood dependence on N . 

In addition, we developed a promising novel way to perform the group integrals approx- 
imately in the limit of large matrix dimension N . Without excessive effort, this method 
improves the traditional Wick pair contraction technique and allows for the calculation of 
corrections in subleading order. It would be of interest to attain a deeper understanding of 
this procedure from a group-theoretical point of view. 
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